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ABSTRACT 

We use numerical TV-body simulations of the Orion Nebula Cluster (ONC) to in- 
vestigate the destruction of protoplanetary disks by close stellar encounters and UV 
radiation from massive stars. The simulations model a cluster of 4000 stars, and we 
consider separately cases in which the disks have fixed radii of 100 AU and 10 AU. 
In the former case, depending on a star's position and orbit in the cluster over 10 
years, UV photoevaporation removes at least 0.01 Mq from its disk, and can remove 
up to 1 Mq . We find no dyn amica l models of the ONC consistent with the suggestion 
of Storzer and Hollenbach ( 1999| ) that the observed distribution and abundance of 



proplyds could be explained by a population of stars on radial orbits which spend 
relatively little time near 9 1 C Ori (the most massive star in the ONC). Instead the 
observations require either massive disks (e.g. a typical initial disk mass of 0.4 Mq) 
or a very recent birth for 1 C Ori. When we consider the photoevaporation of the 
inner 10 AU of disks in the ONC, we find that planet formation would be hardly 
affected. Outside that region, planets would be prevented from forming in about half 
the systems, unless either the initial disk masses were very high (e.g. 0.4 Mq) or they 
formed quickly (in less than ~ 2 Myr) and 8 1 C Ori has only very recently appeared. 

We also present statistics on the distribution of minimum stellar encounter sep- 
arations. This peaks at 1000 AU, with only about 4 per cent of stars having had an 
encounter closer than 100 AU at the cluster's present age, and less than 10 per cent 
after 10 7 years. We conclude that stellar encounters are unlikely to play a significant 
role in destroying protoplanetary disks. In the absence of any disruption mechanism 
other than those considered here, we would thus predict planetary systems like our 
own to be common amongst stars forming in ONC-like environments. 

Also, although almost all stars will have experienced an encounter at the radius 
of the Oort cloud in our own system, this only places a firm constraint on the possible 
birthplace of the Sun if the Oort cloud formed in situ, rather than through the secular 
ejection of matter from the planetary zone. 

Key words: planetary systems - Solar system: formation - open clusters and as- 
sociations: individual: Orion Nebula Cluster - celestial mechanics, stellar dynamics - 
accretion, accretion disks 



1 INTRODUCTION 

The formation of planets from the gas and dust around a 
young star i s thought to take sever al million years (e.g. Lis- 
sauer (1993), Pollack et al. (1996)). During this time, the 



star is unlikely to move far from the region where it itself was 
formed, and its early evolution will continue to be influenced 
by the conditions there. Such regions may vary considerably, 
but the dominant environment for Galactic star formation is 
one containing thousands of young stars in a cluster of high 
density (typically 10 4 stars pc -3 ), and in which the spec- 
trum of stellar masses extends to bright O and B type stars 



- from which they are usually known as OB association s 
dMiller fc Scalo 1978|; |Clarke, Bonnell fe Hillenbrand 2000[). 



The presence of these many other stars may have serious im- 
plications for the nascent planetary system, as close stellar 
encounters can tidally disrupt the protoplanetary disk, and 
UV radiation from the massive stars can destroy it through 
photoevaporation. 

The Orion Nebula Cluster (ONC), a nearby star- 
forming region only a few million years old, is just such an 
environment. Containing over 4000 stars within a volume 
five parsecs across, it nestles at the edge of the giant molec- 
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ular cloud in Orion, about 470 pc away, and its spectacular 
appearance has made it a familiar astronomical image. At 
its centre are the four bright Trapezium stars, of which the 
most massive, 1 C Ori, is by far the dominant source of UV 
radiation in the cluster. 

There is evidence to suggest that the majority of stars 
in the ONC have circumstellar disks. Over 40 have now 
been observed directly with the Hubble Space Telescope, ei- 
ther sillhouetted against the background nebula or embed- 
ded inside a bright ionized envelope of matter (a proplvd) 



( |McCaughrean fe O'Dell 1996fc |Q'Dell fc Wong 1996| ; pally. 



O'Dell Ic McCaughrean 2000). The existence of many more 



can be inferred from an excess in the near infrared contin- 
uum emission from stars in the cluster. Measurements by 



Hillenbrand et al. ( 1998 ) found evidence for disks around 
55-90 per cent of the ir sampl e of ~ 1600 stars, and more 



recently, Lada et al. (in press) have narrowed this fraction 



incorporates two-body regularisation algorithms to handle 
close encounters without smoothing. 



2 COMPOSITION AND DYNAMICS 
OF THE ONC 

Although the ONC is relatively close, the large amounts 
of gas and dust present and the bright emission from the 
Trapezium stars have until recently hidden much of its popu- 
lation, and prevented us from accurately measuring its kine- 
matics. In the past decade, howeve r, deep and high re solu- 
tion studies at optical wavelengths ( Prosser ct al. 1994 ) and 
adaptive opti cal techniqes in the infrared (McCaughrean & 
Stauffer 1994) have greatly improved our understanding of 
the cluster's size and composition 
Hillenbrand 



to 80-85 per cent in their observations. 

In this paper we present the results of JV-bo dy simu- 



Hartmann (1998) detected 3500 stars 
down to mass 0.1 Mq within 2.5 pc of the cl uster centre. A 
subsequent surv ey of the central 0.7 x 0.7 pc 2 ( Hillenbrand & 



20 per cent stars in that region. The overall mass spectrum 
was found to be similar to that of the Galactic field. In this 
work we therefore consider a population of 4000 stars, with 
a mean mass of around 0.5 Mq - though the true population 

is probably higher still 

Jones & Walker ( 1988[ ) measured a (three-dimensional) 
velocity dispersion a = 4.3 ± 0.5kms _1 for about 1000 
stars distributed within 2pc of the centre. Combined with 
a half-mass radius Rh = 1 pc, this gives a crossing time 
T c = 2R h /a w 0.5 Myr. The age of the cl uster is difficult 
to determine, but data from Hillenbrand (1997) and Hil- 



lations of the ONC using Aarseth's NBODY6 code ( |Aarsetb| |Carpcnter 200q ) down to mass 0.02 M found an additional 
200C ) , i ,nd their implications for disk destruction by photoe- 
vaporation and stellar encounters. Such an approach is nec- 
essary to properly calculate the effects of photoevaporation, 
because the mass loss rate from a disk due to an incident 
UV flux is determined by its distance from the flux source. 
The mod el we use for this, due to Jo hnsto ne, Hollenbach 
& Bally ( [1998ft , Storzer & Hollenbac h (|l999|) and references 
in Hollenbach, Yorke & Johnstone ( |2000| ), is in agre ement 
with measurements made by Henney & O'Dell ( 199S| ), who 
found mass loss rates of ~ 4 x 10~ 7 Mq yr -1 from proplyds 
close to 6 1 C Ori. If these objects have spent the whole age 
of the ONC in this environment, their initial masses must 
have been greater than 0.8 Mq, and we might expect to see 
some disks today with a non-negligible fraction of this mass 
present. However mm-wavelength observations of dis ks in 
the ONC indicate ( |Lada et al. 1996| ; |Bally et al. 199S| ) disk 
masses no greater than 0.02 Mq.Q The paradox might be 
resolved if 6 C Ori was born only very recently - at most 
5 x 10 4 yr ago - or if the dynamics of the ONC were such 
that the proplyds seen close to 8 1 C Ori today have spent 
most of their lives elsewhere in t he clu ster, as has been sug- 
gested by Storzer & Hollenbach ( 1999| ). We investigate this 
possibility in simulations in which the cluster undergoes a 
collapse from cold initial conditions. 

The relevance of cluster dynamics for the destruction 
of disks by stellar encounters is even more readily appar- 
ent. A full dynamical simulation is necessary to improve on 
analytic approximations in which stars remain in the same 
density environme nt thr oughout the life of the cluster (e.g. 
Clarke & Pringle (1991)). Until recently however, TV-body 
codes suppressed strongly gravitationally focused encoun- 
ters by smoothing the gravitational field on small scales. In 
the ONC, typical encounters are gravitationally focused at 
impact parameters less than 160 AU, so a proper evalua- 
tion of systems on the scale of interest for planet formation 
could not be made with such codes. In this work we use 
NBODY6, which represents the state of the art in dynamical 
simulation on a commercial hardware platform, and which 



lenbrand & Carpenter (200C) suggest that perhaps 85 per 
cent of its stars are less than 2 Myr old (with a mean age of 
~ 0.8 Myr). The cluster thus seems to be dynamically young 
- a few crossing times old. 

Using these parameters, a rough calculation of the virial 
ratio (kinetic energy / potential energy) of the cluster gives 
a value of about 1.5. Some authors have concluded from this 
that either the ONC must be unbound and expanding, or 
there must be a significant amount of gas or stars - p erhaps 
low-mass binar y companions - p r esent and not seen ( Jones 



& Walker 1985; Tian et al. 1996; Hillenbrand & Hartmann 



1995 ; Kroupa 2000). However, reasonable errors in the obser- 



vational parameters can easily account for an error of over 
50 per cent in this calculation. The density profile of the 
cluster is very much that of a r elaxe d system, to the extent 
that Hillenbrand & Hartmann (1998) were able to fit a King 
model - characteristic of globular clusters many relaxation 
times old - to their data. 

In addition, we note that even if the ONC is unbound 
today, it is unlikely to have been so over much of its past 
evolution. A highly unbound cluster will always appear only 
about a crossing time old, since it expands at some velocity 
v ~ a, and its radius R ~ vt, hence T c as Rja ~ t. Also, 
reasonable initial velocity distributions tend to evolve to a 
density distribution which is much too flat when compared 
to that of the ONC (Scally, in preparation) .£] 



* There are, however, many uncert ainties involved in these esti- 
mates. See Henney & O'Dell (1999). 



t Note that the asymptotic density distribution of a highly un- 
bound cluster is simply mapped from its initial velocity distribu- 
tion, since ultimately r sa vt for each star. Thus an initial velocity 
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3 PROPLYDS AND PHOTOEVAPORATION 

About 150 young stellar objects in the ONC have been ob- 
served with the central star and disk surrounded by a bright 
extended ionisation front, and a tail streaming in the direc- 
tion away from a nearby massive star - usually 6 C Ori. 
The se obj ects are known as proplyds ( after O'Dell, Wen fc 
Hu dl993b), and n recent observations (pally, O'Dell & Mo 



Caughrcan 2000) about 80 per cent of the objects within 



0.14 pc of the Trapezium stars are of this type. The frac- 
tion seems to rise as one looks closer to the cluster centre 



( O'Dell fc Wong 1996 ), though this may be merely a selec- 
tion effect, since their brightness drops in proprtion to the 
square of their distance from 8 1 C Ori, while the background 
nebula fades less rapidly (O'Dell, personal communication). 

A model for these obje cts ha s been developed by John- 
stone, Hollenbach & Bally fll99Sj ) and Storzer & Hollenbach 
( [1999 ) (see Hollenbach, Yorke & Johnstone j2000| ) for a re- 
view), in which their structure and mass loss rate at a given 
distance from a massive star is determined by the relative 
strengths of the star's far ultraviolet (hv < 13.6 eV, here- 
after FUV) and extreme ultraviolet (hv > 13.6 eV, here- 
after EUV) fluxes at that point. Specifically, they propose 



that thjjie aie Iwu itigimus in which pholouvapoialion opei- 
ates. In 1 the inner one, FUV photons ominate the mass loss 

by heating the circumstellar disk and causing a neutral flow 
out to an ionisation front, where it meets the EUV field, 
as seen in the proplyds. Within this regime the model pre- 
dicts a mass loss rate roughly independent of the distance d 
from the massive star, since the only criterion for it to apply 
is that the FUV flux is sufficiently strong to heat the disk 
matter above its escape velocity. In the EUV regime, there 
is no neutral flow and the mass loss rate depends directly on 
the EUV flux - and hence on d. Since all proplyds exhibit a 
stand-off distance between the disk and the ionisation front, 
this model will produce them only in an FUV-dominated 
region. 

Using the same phy sical p arameters as were assumed 
by Sto rzer fc Hollenbach (1999), Hollenbach, Yorke & John- 



stone ( 2000 ) give the following expressions for the mass lost 
by a disk in the FUV and EUV dominated regions: 



M F 



Mf 



2 x lO~ 9 r d M yr~ 1 



8xlO- 12 r^YjM yr- 



(1) 
(2) 



where r<j is the disk radius in AU, 5>i is the ionising (EUV) 
photon luminosity of the massive star in units of 10 49 s _1 , 
d is its distance in pc, and we assume a column density of 
5 x 10 21 cm~ 2 from the ionisation front to the disk inside a 
proplyd. 

For 6 1 C Ori, of type 06 ( |Hillenbrand 199?| ), $i = 2.6, 
and the boundary between the EUV and F UV dominated 



regions stand s about 0.3 pc from the star (storzer fc Hol- 



lenbach 1999). The great majority of proplyds in the ONC 
have indeed been observed within this distance, and it may 
be that those fe w out side (~ 10 per cent in the survey of 
O'Dell & Wong (1996)) involve a wind production mecha- 
nism other than FUV heating. The other bright stars in the 



cluster are all less massive and have significantly less UV 
output, with $i ~ 0.2 for the second most massive (type 
09), and $i ^ 0.05 for the third (type B0). Their FUV- 
dominated regions will be correspondingly smaller (~ 0.1 
pc and 0.05 pc respectively), as will their photoevaporative 
effects. 

The effect on planet formation is complicated by the 
fact that there is a lower limit on the size to which disks can 
be reduced by the photoevaporation process. Ultimately this 
is determined by the mass of the star, since for matter to 
escape from the disk at any point it must be heated above 
the local escape velocity, corresponding to the limitjj] 

r • « ^ (3) 

where M* is the stellar mass and c s is the sound speed in the 
heated flow. For a 0.3 Mq star (corresponding roughly to the 
mode of the mass distribution in the ONC), the minimum 
radius in the FUV-dominated region will be about 15 AU, 
while in the EUV-dominated region evaporation can con- 
tinue down to 1-2 AU - the difference being a consequence 
of the differing sound speeds in EUV an d FUV heated flows 
(lOkms -1 and 3km s -1 resepctively) (Johnstone, Hollen- 



bach & Bally 1995). Even in an FUV-dominated region, 



once the disk becomes too small for FUV radiation to drive 
a neutral flow, EUV-dominated mass loss will set in for that 
system. 

Within the minimum radius, planets may well be able 
to form without hindrance, and for most systems we may 
expect photoevaporation to have very little effect on planets 
forming at about 5-10 AU from their central star (and none 
whatsoever on planets in the inner 1 AU). But this central 
zone is exactly the region we are most interested in for planet 
formation: the planets in the Solar System are found there, 
as are many of the close gas giants that have been discovered 
recently around nearby stars. 



4 SIMULATIONS 

We construct a dynamical mod el of the ONC , implemented 
using Aarseth's NBODY6 code ( Aarseth 200C ), consisting of 
4000 stars, and starting from a density distribution going 
as r~ 2 . Except where stated otherwise (in the discussion of 
Figure ||) , the results shown refer to a cluster in virial equi- 
librium with a half-mass radius of ~ 1 pc, and in all cases 
the initial conditions are chosen to match the appearance of 
the ONC after an evolution time of 2-3 Myr (a few crossing 
times). T he m ass function used is that of Kroupa, Tout & 
Gilmore ( 1995| ), and the three most massive stars are as- 
signed the UV flux parameters specified in the previous sec- 
tion ($i = 2.6, 0.2 and 0.05). Of these, the most massive is 
placed at the cluster centre, while the other two are given 
random initial locations. Several random realisations of this 
setup were generated and run, but no statistically significant 
variations were found, and the results presented here for one 
particular model are characteristic of all those generated. 



distribution uniform in v (say) produces a density distribution 
uniform in r. 



■F The additional factor of ~ 0.5 is due to the fact that material 
can actually e scape from inside GM./c^ (Johnstone, Hollenbach 
fe Bally 1998|. 
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To model photoevaporation, we run two series of simu- 
lations. In the first we are interested in the distribution of 
proplyds (in the FUV-dominated region) and for every star 
in the cluster, so we keep track of three things during the 
course of the simulation: 

• the cumulative time it spends in the FUV-dominated 
region of any UV source. 

• the time-integrated value of yj $i/d 2 for all three UV 
sources whenever it is outside an FUV-dominated region. 

• the closest approach r c i osc it makes to any other star. 

We then use the first two of these in (1) and (2) to evaluate 
the total mass lost by each star due to photoevaporation 
during its life in the cluster by assuming a fixed disk radius 
for all the systems throughout - which we take to be 100 
AU. 

The assumption that the disk radii remain fixed 
throughout is a simplification, since in reality it is likely that 
the disks will decrease in size as they lose mass through evap- 
oration. This would in turn lead to a reduction in the mass 
loss rate, which scales with the disk radius rd in both the 
EUV and FUV regim es, as shown in (1) and (2). Johnstone, 
Hollenbach & Bally ( 199S ) calculate the implied variation 



with time t (at a fixed distance from the ionising star) as 
rd oc t^ 1 (EUV) and rd oc t~ 2 (FUV) for a disk whose sur- 
face density goes as r -3 ' 2 , and assuming that the disk is 
unable to replenish the outer regions from which material 
evaporates. 

Direct measurements of disk radii are difficult to make, 
and to a certain extent depend on the wavelength one ob- 
serves at. However, estimates for various proplyds in the 
ONC ( jjohnstone, Hollenbach fc Bally 1996] ) vary from 20 to 
80 AU, and are in agreement with the model's predictions 
taking into account the size of the ionisatio n front stand-off 



and th e dista nce to 9 C Ori in each case (Storzer&Hol 
lenbach[l99|). Given this, and the possibility that the disks 



might have been much larger in the past if they have already 
suffered a long period of photoevaporation, we take a fixed 
value of 100 AU for these objects as a simplifying assump- 
tion which will underestimate the total mass evaporated to 
date. It would of course be possible to repeat our calculations 
while also keeping track of the disk size for each system, but 
this would necessitate making an assumption about their 
initial size - as well as incorporating the assumption that 
there is no viscous replenishment of the outer disk. This 
would seem to be unwarranted in the light of our finding 
(later in this section) that the mass loss assumption already 
implies unacceptably large initial disk massses for the pro- 
plyds. 

In the second series of simulations we are interested in 
mass loss from the planet-forming region of the disk, i.e. the 
inner ~ 10 AU. We assume that there is no FUV-driven 
mass loss from this region - which is true for any star more 
massive than ~ 0.2 Mq (more than 60 per cent of cluster 
stars) - so that for each star we only need to keep track 
of the time-integrated value of -J <&i/d 2 for all three UV 
sources. This is converted to a mass loss in (2) by setting 
rd = 10 AU. Once again the fixed disk radius assumption is a 
simplification, which in this case will overestimate the mass 
loss from the disk region we are interested in. (We assume 
that the existence of disk material at radii greater than 10 
AU does not increase the mass loss rate within that radius.) 



4000 stars 
2.89 Myr 
100 AU disks 



0.01 0.10 1.00 

Moss loss (M Q ) 

Figure 1. Histogram of mass loss due to photoevaporation after 
2.89 Myr in the case where all disks have radius 100 AU through- 
out. 



4.1 Photoevaporation 

Figure ^ shows a histogram of mass loss due to 
photoevaporation after 2.89 Myr (roughly the ONC's cur- 
rent age) in the case where the disks are 100 AU in radius. 
We see that almost all systems have lost more than 0.01 
Mq, The rightmost bin contains systems which have spent 
their whole time in 6 1 C Ori's FUV domimated region. 

Figure |^ shows, at various times during the life of 8 1 C 
Ori, the percentage of stars in the central projected 0.15 
pc of the cluster which are proplyds, as a function of the 
initial disk mass (assuming all disks are equal initially). A 
proplyd in this context is a star in an FUV-dominated region 
with some circumstellar disk matter remaining (and recall 
that the radius of 8 1 C Ori's FUV-dominated region is 0.3 
pc). At the start of the simulation (or immediately after the 
massive stars formed if that were different), all stars in the 
central 0.3 pc (in 3D) would be proplyds, corresponding to 
<; 80 per cent of the stars within 0.15 pc in projection, given 
the ONC's density distribution. But low mass disks are very 
quickly destroyed, and to match the observed distribution 
of proplyds at the ONC's present age requires a high initial 
disk mass (0.4-0.6 Mq in Figure El) for these stars. (Note 
that these initial disc masses would have to be even higher 
if we were to include the effect of evaporation reducing the 
disk radii and mass loss rates.) 

Figure ^ is a similar plot showing data from a model 
starting in an initially cold configuration, with potential en- 
ergy much greater than its kinetic energy. Such a system un- 
dergoes an initial collapse, after which it virialises at about 
half its initial size (subject to some contnuing small oscil- 
lation in radius). The period of maximum density in the 
collapse is called the crunch, and models which have just 
passed through it can look similar to the ONC. The initial 
density profile is important however, and models starting 
from an r~ 2 profile provide a better match to the ONC at 
its present age than those starting from uniformity. In the 
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Figure 2. Variation with initial disk mass of the percentage of 
stars inside 0.15 pc (in projection) which are proplyds (i.e. in 
an FUV-dominated region with some circumstellar disk matter 
remaining) after 0.96, 1.9, 2.9, 3.9, and 4.8 Myr. All disks are 
assumed equal in mass initially and have fixed radii of 100 AU. 
At 4.8 Myr the percentage has dropped dramatically because we 
assume 9 1 G Ori leaves the main sequence at 4.6 Myr, after which 
only one or two less massive (but longer lived) stars remain to 
ionize any disks. 



latter case, the core tends to be underdense, to the extent 
that only about half the stars appearing in the central 0.15 
pc in projection are within 0.3 pc in 3D. Such a model can 
therefore never match the observations, in which 80 per cent 
of stars in the centre are proplyds. The plot shown here is for 
a cluster after 2.8 Myr with an intial r~ 2 density profile. It 
seems that the dynamics have made little difference to how 
the fraction of proplyds varies with initial disk mass, when 
compared to the case of virial equilibrium. Stars which fall 
into 6 1 C Ori's FUV-dominated region during the collapse 
tend to stay there, and so for low-mass disks, the population 
of depleted disks builds up even as new proplyds appear. 

Note also that as time progresses, for any reasonable 
cluster dynamics, the distribution of proplyds around 6 1 C 
Ori should deplete from the inside out, as these stars are 
more likely to have spent more time in the FUV-dominated 
region. This assumes, of course, that there is no other signif- 
icant mechanism for creating proplyds, and that accretion 
does not replenish the disks. 

Figure ^ shows histograms of mass loss due to 
photoevaporation from 10 AU disks after 2.89 Myr (a) and 
12.5 Myr (b), by which time the massive stars will have left 
the main sequence and photoevaporation will have finished. 
At this point, 95 per cent of systems have lost less than the 
mass o f the minimum Solar n ebula (~ 0.013 Mq (Hayashi 



Nakaza va fc Nakagawa 1985 )) (and recall that this is an 
overestimate of the mass lost from the 1-10 AU region of 
their disks). This is consisitent with the finding of Hillen- 
brand et al. (1998) that in the centre of the ONC there is no 



Figure 3. As Figure |j but for an initally cold cluster after 2.1 
Myr. 



decrease in the incidence of disks revealed by near infrared 
excess. 

4.2 Encounters 

Figure [| shows histograms of r c i osc for all the stars in the 
cluster at the ONC's present age and at 12.5 Myr. In this 
particular run we find that only about 3 per cent of stars 
have had an encounter 100 AU or closer at the present age, 
rising to 6 per cent after 12.5 Myr. Typical values for these 
fractions, averaged over a number of different runs, are ~ 4 
per cent and 8 per cent respectively. 

Figure ^ shows all 4000 stars in the cluster on a plot 
of r c i osc against radial position, both in projection and in 
three dimensions, at the ONC's present age. On the plot in 
projection, we see no trend with i? pr oj outside 1 pc from the 
cluster centre, with almost all stars having a value of r c i osc 
between 10 3 and 10 4 AU, even though the stellar density 
drops by more than an order of magnitude over the range 
plotted. 

It is interesting to compare this with what one would 
expect if each star had spent its whole time at the radius 
observed today - i.e. if no orbit in the cluster had any signif- 
icant radial component. To make this comparison, we per- 
form a Monte Carlo data simulation, the details of which are 
presented in the Appendix. In outline, the simulation consid- 
ers stars in radial bins, each bin having a local density which 
determines the rate of encounters within any given distance. 
Encounters are generated and assigned to the stars in the bin 
at random, and we take the closest encounter for each star 
after the required time (2.89 Myr in this case, to compare 
with the dynamical simulation) . This gives a distribution in 
r c ioso for each radial bin; Figure ^| plots its median and its 
2a and 3er limits against the data from NBODY6. 

The Monte Carlo simulation matches the dynamical 
data remarkably well, with the median, in particular, provid- 
ing a good fit to the 'centre' of the distribution of stars. Its 
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Figure 4. Histograms of mass loss due to photoevaporation after 
2.89 Myr (a) and 12.5 Myr (b) in the case where all disks have 
radius 10 AU throughout. 



Figure 5. Histograms of r c i oso for all the stars in the cluster after 
2.89 Myr (a) and 12.5 Myr (b). 



outer limits in r c i osc are narrower, and this may reflect the 
influence of orbital mixing in the dynamical simulation, but 
the effect is minor, and we conclude that simple analytical 
estimates give a good measure of the encounter distribution 
at 2-3 Myr. 

We find, however, that in evolving from 10 s to 10 7 years, 
the closest encounter distribution hardly changes. This is in 
marked contrast to the analytic estimates, which imply that 
as a result of the increase in available time the distribution 
s hould m ove to smaller radii. For example, Bonnell et al. 
( in press ) find that after 10 7 years, all of the systems in the 



core of a cluster like the ONC (where the stellar density is 
10 4 pc -3 ) have had an encounter within 100 AU, whereas in 
our models less than one third of the systems in the core 
have had such a close encounter. The difference between the 



two is mainly due to the fact that our model clusters be- 
come significantly less dense over this timescale. This evo- 
lution is not driven by mass loss or by super-virial initial 
conditions (though these would have a similar effect), but is 
a consequence of the fact that our initial conditions do not 
represent a steady-state solution of the collisionless Boltz- 
mann equation. The model clusters are generated with a 
sharp discontinuity at the cluster edge, where the density 
drops to zero. In principle, the ONC could be the centre 
of a much more extended distribution of stars comprising a 
steady state distribution, as in the King model fit of Hillen- 
brand & Hartmann ( 199c )|§l and in such a case the density 



S Though such a cluster, to properly match the relevant King 
model, would need to extend to about 20 pc from 9 C Ori and 
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Figure 6. (a) Closest encounter distance against projected radial 
position after 2.89 Myr. (b) Closest encounter distance against 
3D radial position after 2.89 Myr. The thick line is the median of 
a distribution of closest encounters produced by a Monte Carlo 
simulation (see Appendix) and the thin lines and long-dashed 
lines are 2n and 3<r limits respectively on this distribution. The 
short-dashed line is the mean stellar separation at R. 



evolution would not occur and the core encounter distribu- 



tion predicted by Bonnell et al. (in press) would be appro- 
priate. However, we point out that for initial conditions that 
are not so carefully constructed, we expect the core density 
to fall and for most of the closest encounters in the ONC to 
have occurred by its present agej^] 



5 DISCUSSION 
5.1 Proplyds 

It is clear from the mass loss rate in (Q) that photoevap- 
oration can remove a large amount of disk material from 
systems in the core of the ONC over the ~ 5 Myr lifetime 
of 8 1 C Ori. Our simulations, which keep track of the stars' 
varying distances from 1 C Ori as they orbit in the clus- 
ter, allow us to quantify the fraction of systems that never 
spend any significant time near 9 1 C Ori, and remain rela- 
tively unscathed by photoevaporation. They therefore shed 
some light on the much discussed problem of why proplyds 
are common in the core of the ONC, when the short inferred 
survival times of disks in this region would imply that most 
stars have already have lost their disks and hence should not 
manifest proplyd activity. 

One solution, proposed by Storzer and Hollenbach 

([l99§), is that the proplyds are merely 'visiting' the core 
region on radial orbits. In this picture, stars light up as pro- 
plyds inside the cluster core, but since they have spent most 
of their lives at much larger radial distances, their disks can 
have survived for much longer. We find no evidence for such 
a population of stars on radial orbits in any dynamically 
plausible model for the ONC. In models that are in virial 
equlibrium, the velocity dispersion is initially isotropic, and 
remains approximately so, and the stars currently in the 
central region have spent most of their lives there. We find 
that in this case the observed high proplyd fraction is con- 
sistent only with initial disk masses in excess of 0.4 Mq 
for somewhat smaller initial disk masses, the proplyd dis- 
tribution develops a ring shaped structure in projection, as 
only disks towards the edge of the FUV-dominated region 
(which spend some time outside it) can still exist. The lack 
of any observational evidence for such a central depletion in 
the proplyd distribution suggests that proplyds are not close 
to the point of exhaustion. 

Perhaps more remarkably, we find that the 'radial orbit' 
solution of Storzer and Hollenbach also fails to be realised 
in the case where the cluster undergoes cold collapse. In 
such models (where the kinetic energy of the cluster is ini- 
tially low, as would be the case if the stars had fragmented 
out of a hydrostatically supported medium) the stars begin 
by falling inwards on radial orbits. After about a free-fall 
time, they achieve a configuration of maximum compact- 
ness (the crunch) and rebound into a state of approximate 
virial equilibrium. At an age of about 2 Myr, the ONC would 
have evolved somewhat past the crunch and its orbits would 
be largely isotropic. In consequence, the predicted proplyd 
fraction is remarkably similar to the virialised case (Figure 
3 compared with Figure 2). 

Since we have shown that orbital dynamics cannot solve 
the proplyd frequency problem, we are forced to invo ke the 
other two s o lutions considered by p revious authors (Bally 



et al. 1998 



Henney fc O'Dell 1999). The first possibility 



is that the measured disk masses are vast underestimates. 



contain about 30,000 stars, which is almost certainly not consis- 
tent with observations of the Orion region. 

» Another dynamical factor not present in the Monte Carlo sim- 
ulation is that for stars in the cluster core, many close encoun- 



ters may occur within small- ./V groupings, rather than being dis- 
tributed over the entire core population. 

II Which would of course constitute an unstable disk fo r star s 
of Solar mass or less (see e.g. Laughlin & Bodenheimer (1994), 
Toomre (1964)). 
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In order to solve the proplyd frequency problem, the initial 
disk mass required (0.4 Mq) would imply that 90 per cent 
of the stars only lose a small fraction of their disk mass to 
photoevaporation over the lifetime of 6 1 C Ori. The second 
possibility (which would be compatible with the low mea- 
sured disk masses) is that 6 1 C Ori has formed only recently 
(implying incidentally that the ONC has been caught at a 
special moment, and that older clusters containing OB stars 
would not be expected to exhibit proplyd activity). 



5.2 Planets 

One might suppose that the high mass loss rates implied 
by (1) should lead to a suppression of planet formation in 
all populous cluster environments containing O stars, unless 
planets form at the same time as the disk itself. Certainly 
this is the case for planets forming outside the inner 10 AU 
of the disk, since as Figure |l] shows, even after 3 Myr most 
systems with 100 AU disks will have lost an amount of disk 
material equal to several times the minimum Solar nebula. 
Planets would be prevented from forming on such wide or- 
bits in about half the systems in the ONC, unless, as dis- 
cussed above, either the initial disk masses were very high 
(e.g. 0.4 Mq) or they formed quickly (in less than ~ 2 Myr) 
and 8 1 C Ori has only very recently appeared. 

However, it is the inner 10 AU where we expect most 
planets to be found (by analogy with our own system) and 
as explained in Section 2, only EUV-dominated mass loss 
can affect this region of the disk. Figure ^ shows that such 
mass loss is low - relatively few disks spend much time very 
close to an O star. The results indicate therefore that planet 
formation in the ONC would be largely unaffected by pho- 
toevaporation. 

Disruption due to stellar encounters is probably also 
unimportant. In a star-disk encounter, matter can be 
stripped from the d isk down to about one third of the en- 
counter separation ( Clarke & Pringle 199S ) , and during the 
lifetime of the cluster only a small minority of stars in our 
simulations will have had encounters close enough to affect 
the planet-forming region of their disks. This is essentially in 
agreement with the results of Bonnell et al. (in press) , who 
find that only in the dense core of the ONC is there likely 
to be any noticeable disruption to young planetary systems. 

Our conclusion, therefore, is that in the absence of any 
disk disruption mechanism other than those considered in 
this paper, we would predict planets in orbits within 10 AU, 
and perhaps planetary systems like our own, to be com- 
mon amon gst st ars forming in ONC-like environments. As 
Armitage ( 2000[ ) has pointed out, in much richer stellar en- 
vironments, where the number of O stars is greater, photoe- 
vaporation may have a m ore signi ficant impact. The obser- 
vations of Gilliland et al. ( frn press ) , who find a complete ab- 
sence of close ('hot') Jupiter-mass companions in the globu- 
lar cluster 47 Tucanae (which contains more than 10 5 stars) , 
may be an example of this. More challenging, perhaps, for 



our understanding of star and planet formation are the re- 
sults of the gravitational lensing survey of Albrow et al. (in 
press), who find that less than a third of lensng stars (typi- 
cally of mass ~ 0.3 Mq) have Jupiter-mass companions with 
orbits in the range 1.5-4 AU. If the majority of these lens- 
ing stars were formed in an environment no richer than the 



ONC then neither photoevaporation nor stellar encounters 
can explain the apparent absence of planets. 



5.2.1 The Oort cloud 

We note finally that all stars in the ONC should have had 
an encounter within a radius comparable to that of the Oort 
cloud in our own system (<; 20,000 AU) , and that if the Oort 
cloud is primordial (e.g. Cameron (1973)), this could be used 
to rule out an origin for the Solar system in an ONC-like 
environment. On the other hand, it is often supposed that 
the Oort cloud was formed from bodies scattered out of the 
planet-forming zone by pla netary and tidal perturbations 



ort 1950 



Fernandez 1985), arriving at large radii only af- 



ter about 10 years (Duncan, Quinn & Tremaine 1987). If 



the ONC is unlikely to survive as a bound cluster for that 
duration, its density would be considerably less by then, 
and an Oort type cloud might be able to survive without 
difficulty in the more dilute environment. Given the uncer- 
tainties in the origin of the Oort cloud, we conclude that its 
existence alone is not sufficient to place any firm constraints 
on the birthplace of the Sun. 
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stars at radius r in the cluster. Taking the closest encounter 
for each star then gives a distribution of minimum encounter 
distances at r, analagous to that obtained from the dynam- 
ical simulation (with which we match the parameters k, a, 
m and T) . Figure |^ plots the median of the distribution and 
its 2a and 3a limits. 



APPENDIX A: MONTE CARLO SIMULATION 
OF ENCOUNTERS 

In an environment with a given local stellar density n and 
one-dimensional velocity dispersion a, the rate of encounters 
(per star) closer than a given distance r e can be calculated 

as 

/ = 4y^n (arl + (Al) 

(e.g. Binney & Tremaine (1987), p. 541) where m is the 
stellar mass, assumed the same for all stars. In a cluster 
where the density n(r) = kr~ 2 for some constant k there 
are 4nk Sr stars with radial positions in [r, r + Sr] , and after 
a time T we expect them to have 

N c (r,5r,r c ) = lQn 3/2 k Sr n{r)T (ar 2 + (A2) 

encounters in total closer than r c . Normalising this to 
N e (r,8r,R) and inverting gives a generating function for 
stellar encounters closer than some maximum distance R: 



r c = \J A 2 + BX - A (A3) 
where X is a uniform random deviate in [0, 1] and 

^=£% B = R 2 + ^ (A4) 

Gm 

For R we take the mean stellar separation at r: 

1/3 



R = — - - (A5) 
ynn(r) J 

since all stars have encounters at this distance or greater at 
any given instant. 

Our procedure is then to generate N e (r,5r,R) encoun- 
ters using (|A3|) and assign them randomly to the AnkSr 
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